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Abstract. An important question in the study of N — 2 supersymmetric string or field 
theories is to compute the jump of the BPS spectrum across walls of marginal stability in the 
space of parameters or vacua. I survey four apparently different answers for this problem, 
two of which are based on the mathematics of generalized Donaldson- Thomas invariants 
(the Kontsevich-Soibelman and the Joyce-Song formulae), while the other two are based 
on the physics of multi-centered black hole solutions (the Coulomb branch and the Higgs 
branch formulae, discovered in joint work with Jan Manschot and Ashoke Sen pQ). Explicit 
computations indicate that these formulae are equivalent, though a combinatorial proof is 
currently lacking. 
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1. Introduction 

In quantum field theories and string theory vacua with extended supersymmetry, it is often 
possible to determine the spectrum of BPS bound states in some weakly coupled region of 
moduli (or parameter) space B. In extrapolating the BPS spectrum to strong coupling, one 
usually faces two issues: i) short BPS multiplets may pair up into long multiplets and leave 
the BPS spectrum and ii) single-particle bound states may decay into the continuum of multi- 
particle states. The first issue can be avoided by considering a suitable index 0(7, t), designed 
such that contributions from long multiplets cancel. 0(7, t) is then a piecewise constant function 
of the charge vector 7 and couplings/moduli t G B. 

The second issue arises at certain loci in moduli space, where the bound state becomes 
unstable towards decay into a re-particle state with charges {cti} such that 7 = J2i=i n a i- ^ n 
four-dimensional field or string theories with J\f = 2 supersymmetry, the mass of a BPS bound 
state M(7, t) is equal to |Z(7,i)|, where the central charge Z is a map from B to Hom(r, C), 
where T is the charge lattice. In particular, Z is linear in its first argument 7. The decay is 
therefore energetically possible only when (and even then, marginally so) the phase of Z(~f, t) 
aligns with the phase of each of the Z(ai,t) , s, so that M(j) = ^i=i n-^( a «)- This alignment 
takes place in a locus of codimension p — 1 in the moduli space B, where p is the dimension of 
the subspace of T spanned by the a^s. The dangerous case is p = 2, where the locus defines a 
codimension one "wall of marginal stability" in B, across which the index f2(7,i) may jump. A 
paradigm of this phenomenon is Seiberg-Witten theory with SU (2) gauge group and no flavors: 
across the curve {ajan G M + } in the u-plane, the BPS spectrum jumps from an infinite number 
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Figure 1. Chamber structure of the u-plane and BPS spectrum in N = 2, D = 4 SYM theory 
with SU(2) gauge group and no flavor. The line Im (a/a£>) = separates the strong and weak 
coupling chambers. The only stable BPS states in the strong coupling chamber are the monopole 
and dyons with charges (q,p) = ±(0, 1), ±(2, —1), in the conventions of [5]. The weak coupling 
spectrum consists of these same states and their images around the monodromy at infinity, plus 
the VF-boson with charge (2,0). 

of states in the weak coupling region to just two states in the strong coupling region - the 
monopole and the dyon (see Fig. [I]) [2 [3]. An important physical question is therefore to 
determine the jump 

Aft(M 7 i + iV 72 ) = 0~(M7i + iV 72 ) - ^ + (M 7 i + iV 72 ) (1-1) 

of the BPS index Qo( 7 ,i) = Tr'(-1) F across the wall of marginal stability 

^(71,72) = {t£B/ arg[Z( 7l ,t)] = arg[Z( 72 ,t)]} , (1.2) 

in terms of the BPS indices on side of the wall, say the fi + ( 7 )'s. Here we denoted by ^^(7) 
the index in the chamber c on the side of the wall where argZ( 7 i,i) ^ argZ( 7 2,t). It will 
be convenient to choose the basis 7 i, 7 2 such that + (M 7 i + ./V72) vanishes whenever MN < 
(the 'root basis' condition [4]), and denote by T = (Z + 7 i + Z + 7 2)\{0} the positive cone in the 
two-dimensional sublattice spanned by 71,72- We further assume that ( 7 i, 7 2) < 0. 

As we shall explain in § |4.2[ in M = 2 supergravity and for suitably large charges, the jump 
in f2(M 7 i + iV 7 2,i) is accounted by the loss or gain of a family n-centered BPS black hole 
solutions with charges cti = Mj 7 i + iVj 7 2 G T, with X^(Mj, iVj) = (AI,N), which exist on the 
side c_ of the wall, and on this side only. Close to the wall, the n-centered configuration is 
loosely bound, with the relative distances between the centers diverging at the wall. As 
a result, its index factorizes into the product of the internal index fi(aj) associated to each 
center, times the index of the degrees of freedom associated to the relative motion of the centers 
(with suitable modifications due to Bose- Fermi statistics when some of the (Mi, Ni) coincide) [Ij. 
In §4.2[ we shall compute this configurational index using localization techniques, and obtain 
AO(M 7 i + iV 7 2) for arbitrary values of (M,N). 

A closely similar problem arises in the mathematics of Donaldson-Thomas invariants of 
coherent sheaves on a compact complex manifold X. These invariants, which we shall again 
denote by f2( 7 , t), are labelled by a class 7 in the K-theory lattice V = K(X), and depend on a 
choice of stability condition <fit '■ T — > S 1 inside a complex family parametrized by t. 0( 7 ,t) is 
defined as the Euler characteristic P] of the moduli space A^( 7 , t) of stable coherent sheaves on X 

1 Here, Tr ' denotes the trace in the Hilbert space associated to (7, t) with the center of motion degrees of freedom 
removed, and (— 1) F denotes the fermionic parity, equal to (— l) 2 '' 3 by the spin-statistics relation, where J3 is the 
angular momentum operator along the z axis. 

2 Rather, the Euler characteristic weighted by Behrend's function, see ^\ [7] for details. 



in the class 7 with respect to (j)f As the stability condition is varied, some of the stable sheaves 
may become unstable and the DT invariant 0(7, t) may jump. This happens on the same walls 
of marginal stability as in ( 1.2 ), where ^(7) plays the role of arg Z{p/, t). The similarity between 
these two wall-crossing problems follows from the fact that stable objects in the derived category 
of coherent sheaves on a Calabi-Yau three-fold X are realized physically by BPS states in type 
IIA string theory compactified on X (see e.g. 019] for reviews). 

In two independent pieces of work, Kontsevich-Soibelman [10J and Joyce-Song [7J have 
determined the variation in terms of the DT invariants on one side of the wall. In 

both works, it was noted that the wall-crossing formula takes a simpler form in terms of the 
'rational DT invariants' ^(7, t), related to the ordinary, integer-valued invariants ^(7, t) by the 
'multi-cover formula' 

n(7,t) = ^0( 7 /m,t)/m 2 , (1.3) 
m|7 

where the sum runs over all integers rn > 1 such that 7/m G T (thus, ^(7) = ^(7) if 7 
is primitive). From a physics point of view, we shall explain in £ ]4.1| that the replacement 
0(7, t) —7- t) effectively converts the Bose-Fermi statistics of the n centers into Boltzmannian 
statistics, thereby allowing us to treat the centers as distinguishable. In more detail, the KS 
and JS formulae express the jump Ad (7) as 

Aii(7. f ) = £ £ fg?n:- i n + «-»')> M 

«>2 { ai ,...a n }ef 11 
7=«iH \-a n 

where the second sum runs over all (unordered) decompositions of the total charge vector 7 
into a sum of n vectors ojj G V. The coefficients g({aii}) (after extracting out a Boltzmann- 
Gibbs symmetry factor Y\ m^- whenever {on} contains m\ copies of f3\, m2 copies of 02, etc) are 
universal functions of the a^s. As we shall see, the function g{{at}) turns out to be equal to 
the index Tr'(— 1) F of the configurational degrees of freedom of an n-centered Boltzmann black 
hole solution. To compute this index, it is convenient to consider the refined (or equivariant) 
index g({oti}, y) = Tr '(— y) 2J ' A , evaluate the latter by localization methods, and set y = 1 at the 
end. In fact, the refined configurational index g({cti},y) for general y enters the wall-crossing 
formula for so-called motivic Donaldson-Thomas invariants, see (2.11). 

In the remainder of this survey, we first give an executive summary of the KS (Q and JS 
(§[3]) wall-crossing formulae, silencing the subtleties involved in defining the Donaldson-Thomas 
invariants by themselves. In ^4] we then derive the combinatorial factors g({cti}) by quantizing 
the phase space of multi-centered BPS black holes and evaluating the index by localization. In 

we give an alternative computation of g{{a{\) relying on Reineke's results for quivers without 
closed loops. We end in ^6] with a discussion of some open problems. The material is mostly 
based on [HE], which the reader should consult for more details. Other important references 

include nanamiiaiigiiaEiiniiisi. 



2. The Kontsevich-Soibelman wall-crossing formula 

We start by reviewing the Kontsevich-Soibelman wall-crossing formula for generalized 
Donaldson-Thomas invariants. For generality, we first state its motivic (aka refined) version, 
and then discuss its classical limit. As an application, we use the KS formula to rederive the 
primitive and semi-primitive wall-crossing formula. Finally, we extract the coefficient g{{ati}) 
appearing in (1.4) for n < 3, for comparison with other wall-crossing formulae discussed later. 



2.1. The motivic KS formula 

The motivic KS formula pertains to 'motivic Donaldson- Thomas invariants' Vt n {li t) attached 
to a Calabi-Yau threefold category X, with a stability condition (fit- Informally, Q n (j,t) is the 



n+hd-ih Betti number of the moduli space M. (7, t) of stable degree-7 objects in the triangulated 
category of coherent sheaves on X, where d = dim At (7, t). We define the Poincare polynomial 



9. 



rcf ( 7 ,t,y) = ^(-y)"ft n ( 7 ,t) 



(2.1) 



net 



which is a finite Laurent polynomial in y, symmetric under y — >■ 1/y (the subscript 'ref stands 
for 'refined', which for our purposes is synonymous with 'motivic' |17j). 

To state the KS formula, we introduce the Lie algebra A spanned by abstract generators e 7 , 
for each 7 G K(X) = T, subject to the commutation rule 



7, e 7'J 



^((7,7')>2/)e 7 + 7 ', 



(2.2) 



where (7, 7') is the integer-valued antisymmetric pairing on T (physically, the Dirac-Schwinger- 
Zwanziger product for electromagnetic charge vectors), and 



n(x,y) 



{-y) x - (-y)~ 
y - i/y 



. sinh(z^x) 
sinh v 



\sxy 



(2.3) 



It is straightforward to check that (2.2) satisfies the Jacobi identity for any y. 

Now, for a given choice of stability condition t and any 7 € T , let U^(t) be the element in 
the group Q = exp(^4) defined by 



u 7 (t) = n e 



y n e 7 



y — i/y 

where E(x) is the quantum dilogarithm function 



(-1)™ +I n„(7,t) 



(2.4) 



E(x) = exp 



(xy) k 



^ k(l -y 2k ) 



(2.5) 



We now restrict to 7 G f , where T is the positive cone in the two-dimensional sublattice of 
T spanned by two primitive vectors 71,72, and assume that the 'root basis' condition stated 
below (1.2) holds. The motivic KS wall-crossing formula \10 \ 1171 [18] states that the following 
two ordered products 



n 



u 



M71+N72 



n 



u 



M7i+Af72 ' 



(2.6) 



M>0,N>0, 
M/Nl 



M>0.N>0, 



where the products are ordered with decreasing (resp., increasing) values of M/N G [0, +00] 
(such that the argument of Z(a) decreases from left to right on either side). Here, denote 
the group element C/ 7 (i) when t lies on the side c± of the wall W(7i,72). Thus, assuming 
that the fi+(7)' s 

are known for all 7 G T, the O n (7)'s can be computed by re-ordering the 



product on the l.h.s. of ( |2.6[ ) in the opposite order, using the commutation rule (2.2) and the 
Baker-Campbell-Hausdorff (BCH) formula, and reading off the exponents ^"(7). 

This procedure is vastly simplified by expressing the products in (2.6) in terms of the 'rational 
motivic invariants' 







ref 



(.1,t,y) = J2 



(y-y x ) 



m\-y 



m(y r 



y~ 



(2.7) 



The relation between f2 re f (7, y) and f2 re f (7, y) is easily inverted by means of the Mobius formula, 



m|7 



(y-y l ) 

m(y m — y~ m ) 



^ref (7 /ni>,t,y r 



(2i 



where fi(d) is the Mobius function (1 if d is a product of an even number of distinct primes, — 1 
if d is a product of an odd number of primes, or otherwise). Using the fact that the generators 
eg-y commute for all £ G Z, we may rewrite ( |2.9| ) as a product of factors labelled by coprime 
(M,N), „ 

(2-9) 



n 



V, 



M7i+iV72 



n 



M71+7V72 ' 



M>0,iV>0>0, 
gcd(M,N)=l,M/Nl 



M>0,N>0>0, 
gcd(M,A r )=l,M/iVt 



where 



(2.10) 



£>i 



\JV=1 



Due to the fact that the algebra (2.2) is graded by T, and that every factor of e a in the logarithm 
of is multiplied by a factor of $7 rc f(a,y) for the same vector a, it is clear that the result of 
the re-ordering procedure outlined below ( |2.6| ) will produce a wall-crossing formula of the form 

5ref({"i},y) n" 



^cf(7,y)-^ + ef(7,y) = E £ 



"> 2 {»!,.. .,a n }ef" 
7=«lH h«„ 



Y\m k l 



II J=1 ^efKy)> (2-n) 



where the second sum runs over unordered sets of n charge vectors ai = M^i + such 
that Yli=i n a i = 7 = -^71 + ^72- By the 'root basis' property, only a finite number of terms 
appear on the right-hand side. In §2.2| below, we shall compute the universal combinatorial 
factors <7 re f ({««}, y) in selected cases. 

Before doing so however, we discuss the classical (or numerical) KS wall-crossing formula, 
which arises from the motivic formula (2.9) in the limit y — ¥ 1. In this limit, the Poincare 
polynomial (2.1) reduces to the Euler-Behrend characteristic of the moduli space At(7,t) 



ref ( 7) t,y) ^ x[M(-y,t)]€Z, 
while fj r ef(7, y) reduces to the 'rational DT invariant' 



(2.12) 



m|7 



(2.13) 



Moreover, the commutation rule (2.2) in the Lie algebra A has a smooth limit 



^ ey] = (_i)<7,7'> (7) y), 



-7+7 ; 



(2.14) 



and so does the operator U 7 in (2.10). As a result, the combinatorial coefficients y re f ({««}, y) 
have a smooth limit as y — > 1, and the wall-crossing formula for the rational DT invariants is 
given by the limit of (2.11) as y — > 1, i.e. Eq. (|1.4[) with 



9{{ui}) = limy ref ({a i },y) 



(2.15) 



2.2. Primitive and semi-primitive wall-crossing 

Despite the fact that either side of the KS wall-crossing formula (2.9) involves an infinite number 
of factors, the procedure of re-ordering the product involves only a finite number of operations, 
for the following reason (already hinted at below (2.11)): for any M,N > 0, 

Im,n = { K • e m7l+n72 } (2-16) 

m>M and/or n>N 

is a two-sided ideal in A, and the quotient Am,N = A/Im.n is a finite dimensional algebra. For 
the purpose of computing Af2(M7i + -/V72), it is sufficient to project the relation (2.9) to Am,N 
and use the truncation of the BCH formula at order min(M, N). E.g. to compute ACl(^y% + 72), 
it suffices to re-order the l.h.s of the identity in Ai t i 

exp(0 + (7i)e 71 ) exp(ft + (7i + 7 2 )e 7l+72 ) exp(0 + (7 2 )e 72 ) 

= exp(fT(7 2 )e 72 )exp(fr(7i + 7 2 )e 7l+72 ) exp(fT(7i)e 7l ) (2.17) 

using the truncated BCH formula e x e Y = e x+Y+ ^ X)Y \ and match the result to the r.h.s. In 
this way, we find that the motivic invariants O xe f (71, y), O re f (72, y) are constant across the wall, 
while 

A0 re f(7l +72,y) = K((71 ! 72),y)^ref(7l>y) fi ref(72,y) • (2.18) 

This relation (or its obvious classical limit at y = 1) is known as the 'primitive wall crossing 
formula' [HI [H2 Q2] . From (2.18), setting oi = 72 and 02 = 71 so that (0:1,02) > 0, we read off 
the combinatorial factor 

5ref («i,« 2 ,y) = -<(aua 2 ),y) = (-!)<— >+* sinh H^)) . ( 2 .19) 

smhi/ 

Up to a sign, this is recognized as the character Tr (—y) 2j3 of a representation of SU(2) with 
spin j = 5(1(01,02)) - 1). 

With some more work, one can easily extract the combinatorial coefficients <? re f ({oj}, y) for 
n > 2. E.g, for n = 3 and 01,02,03 three distinct (non necessarily primitive) elements of V 
ordered such that Ojj = (o,, Oj) > for i < j, we find 

fcf (01, 02, 03, y) = (-l)-^3 +Q23 sinh(,o 12 )sinhKoi3 + o 2 3)) (2 2Q) 

sinh v 

when 012 > 023, or 

^(01,02,03,,) = (-1)— +« 23 sinh(,023)sinhKoi2 + o 13 )) 

Sinn v 

when 012 < 023. The result for n = 4 can be found in [I]. 

While the amount of work necessary to extract the combinatorial factors quickly grows with 
(M,N), for fixed (small) M it is possible to compute all the jumps for 7 — > M71 + ./V72 at 
once using the Hadamard lemma log(e x Y e~ x ) = ^2 n>0 Ad 1 ^ ■ Y)/nl, where Adx ■ Y = [X, Y]. 
Defining 

00 

Z ± {M, q,y) = Y J nf ci (M, N, y) q N , fi± f (M, N, y) = (M 7l + N l2 ,y) , (2.21) 

N=0 

we find, for M = 1 [20|, fT] . the 'semi-primitive' wall crossing formula [T3] 



Z (l,q,y) = Z + (l,q,y) Z heLlo (j!, q,y) 



(2.22) 



where 

Zhalo(7i>Q)2/) = exp ^^«(^(7l>72)>2/)firef(^72,2/)^J • (2.23) 

The reason for the subscript 'halo' will become apparent in §4.1| Re-expressed in terms of the 
integer motivic invariants, this can be written as an infinite product [18J 



Z halo (^, q ,y) = [J i 1 ~ (-l) fe l^l 9 V +2i - 1 - fe|wl ) (_1)nn " (fe72) (2-24) 

k>l,neZ 
l<J<fc|7ia| 



Generalizations of (2.22) for M = 2,3 can be found in [T]. 



Exact wall crossing 

Finally, we discuss some examples where the re-ordering of the product in (2.9) can be performed 
in the full untruncated algebra A. Suppose that in the chamber c+, the only non- vanishing DT 
invariants are Q + (-~fi) and r2 + (72). If 712 = — 1, the result of the re-ordering givea^] 

U 12 V lx = L7 7l U 11+12 lV 72 , 712 = -1 , (2.25) 

which follows from the pentagonal identity for the quantum dilogarithm function. If instead 
712 = — 2, one arrives at [TO] 

^(2,-1) • ^(0,1) = ^(0,1) - ^(2,1) • ^(4,1) • • • £^(2,0) • • • ^(3,-1) ' ^(2,-1)^(1,-1) > (2.26) 

where we denoted 72 = (0, 1), 71 = (2, —1) to match the usual basis of electromagnetic charges 
in Seiberg-Witten theory with G = SU(2) and no flavors [5]. As first noted by Denef, Eq. 



(2.26) then embodies the BPS spectrum of this gauge theory on the two sides of the curve of 
marginal stability Im {a/ajj) = 0, see FigJTJ Analogues of ( |2.26 ) for SU{2) gauge theories with 



< Nf < 4 flavors can be found in |15t I18[7 More general identities of this type can be derived 
using Y-systems and cluster algebra techniques, see e.g. 



3. The Joyce-Song wall-crossing formula 

In this section, we briefly review the Joyce-Song wall-crossing formula, which was derived in the 
context of the Abelian category of coherent sheaves on a Calabi-Yau three- fold X [7] . Unlike the 
KS formula, the JS formula only applies to the jump of the classical (or numerical) Donaldson- 
Thomas invariants. Moreover, it gives a fully explicit formula for the combinatorial factors 



g({ai}) appearing in (1.4). The price to pay is that the JS formula is computationally less 
efficient, as it involves sums over many terms with large denominators and large cancellations. 

To state the JS formula, we first introduce S, U and C factors, which are functions of an 
ordered list of charge vectors aj € T,i = 1 . . . n: 

• We define S(ai, . . . ,a n ) G {0,±1} as follows. If n = 1, set S(u\) = 1. If n > 1 and, for 
every i = 1 . . . n — 1, either 

(a) (ai, aj+i) < and (a± H \-a>i,a i+ i-\ h a n ) < 0, or 

(b) (a i ,a i+ i}>0 and (a± H \-ai,a i+ i-\ h a n ) > , (3.1) 

let S(ai, . . . , a n ) = (— l) r , where r is the number of times option (a) is realized; otherwise, 
S(a\, . . . , a n ) = 0. 



3 In ( |2.25[ ) and ( 2.26[ ), Q, = 1 for each factor, except for the factor (7(2, 0) hi the middle of ( 2.26 1, for which SI — —2. 



To define the U factor (not to be confused with the operator U of the previous 
section !), consider all ordered partitions of the n vectors on into 1 < m < n packets 
{ota.j-i+1, ■ ■ ■ , ce aj }, j = 1 • • • tti, with = ao < a\ < • • • < a m = n, such that all vectors 
within each packet are collinear. Let 



Pj = a aj _ 1+ i H h a a . , j = 1 . . . m 



(3.2) 



be the sum of the charge vectors in each packet. Next, consider all ordered partitions of the 
m vectors f3j into 1 < I < m packets {/3fc fe _ 1 _)_i, • • • , /?&,.}, with = 60 < &i <•••<&/= m, 
k = 1 . . . I, such that the total charge vectors 5k = f3b k _ 1 +i + • • • + f3b k , k = 1 . . . I are all 
collinear. Define the [/-factor as the sum 



U(ai, ...,a n ) = j 

1 

m 1 

n 



d-i 



n 



(3.3) 



over all partitions of on and j3j satisfying the conditions above. If none of the on are 
parallel, S = U. Contributions with I > 1 arise only when can be split into two (or 
more) packets with the same total charge, e.g. 



1 



^[71,72,71,72] = S[7i, 72,7i> 72] - o' S '[7i,72] 



1 - "(-I) 2 = " 

2 V ; 2 



Finally, departing from the notations in [7], define the C factor by 

trees edges(ij) 



(3.4) 



(3.5) 



where the sum runs over all labelled trees with n vertices labelled {1, . . . ,n}, with edges 
oriented from i to j if i < j. There are n n ~ 2 labelled trees with n-vertices, which can be 
labelled by their Priifer code, an arbitrary sequence of n — 2 numbers in {1, . . . n}. 

With these definitions, the result of [7] can be stated as an explicit formula for the 
combinatorial factors g({a{}) appearing in (1.4): 

9{{<*i}) = (-ir~ 1+ ^<>^> Yl C KW> • • • a <n)) U (a CT( i), . . • a a(n) ) . (3.6) 

As an illustration, we now use the JS formula to derive the combinatorial factors g({cti}) 
for n = 2,3. For n = 2 (assuming as before that 712 < 0), the U,S,C factors are given in the 
following table 



(7(12) 


S 


U 


c 


12 


a 


-1 


712 


21 


b 


1 


-712 



The JS formula (3.6) then leads to 

5(71,72) = (-1) 712 712^(71)^(72) 



(3.7) 



(3i 



in agreement with the classical limit of (2.19). 



For n = 3, assuming as before that ai,a2,as are three distinct (non necessarily primitive) 
elements of T ordered such that ay = (a:j,ay) > for i < j and moreover that ayi > a 2 3, we 
find that the S, U, C factors are given by 



a(123) 


S 


u 


c 


123 


bb 


1 


ai2«13 + "13023 + «12a23 


132 


b- 





ai2«13 - ai3"23 - ai2«23 


213 


ab 


-1 


-ai2«23 + «13"23 - ai2"l3 


231 


-a 





012^13 - "13023 - ai2«23 


312 


ab 


-1 


ai3"23 - "12023 - ai3«12 


321 


clcl 


1 


ai3«23 + ai20l3 + «12tt23 



(3.9) 



The JS formula (3.6) then leads to 

g({ ai , a 2 , a 3 }) = (-lj^+aas+aia ai2 ( ai3 + 



(3.10) 



in agreement with the classical limit of (2.20). The computation for n = 4 can be found in [I], 
and matches the result from the KS formula. 

4. Wall-crossing from multi-centered quantum black holes 



In this section, we give a new physical derivation of the wall-crossing formula ( 1.4 ) (in particular, 
a new formula for the combinatorial factors g({cti, y})) based on the quantum mechanics of multi- 
centered black hole configuration. Before starting, it should be noted that the KS wall-crossing 
formula (and to a lesser extent, the JS formula) has already been derived or interpreted in various 
physical settings [191 13 EH ES EH E3 EH 12 EH] . Our derivation is arguably more elementary, 
as it relies only on the supersymmetric quantum mechanics of point particles interacting by 
Coulomb and Lorentz-type forces. The down-side is that it does not make the algebra A manifest 
and, admittedly, relies on some plausible but not rigorously proven assumptions. 

4-1- From Bose-Fermi to Boltzmann statistics 



To motivate our approach, let us return to the semi-primitive wall-crossing formula (2.22) and 



rewrite (2.22) as 



for simplicity, concentrate on the classical limit y = 1. Substituting (2.13) in (2.23), one may 

(4.1) 



N 



E.g. for 7 i— > 71 + 272, we find 



n (1 - (-i) fc71 

fc>0 



fcl 712 l n+(fc 72 ) 



ao(i, 2) =2 7l2 n+(i, 0) n + (o, 2) + (-1) 712 7i2 n+(i, i)n+(o, 1) 

"1 



+ ^ + (l,0) 



- 7l2 O + (0,l) ( 7 12^ + (0,1) + 1) 



(4.2) 



The two contributions on the first line can be interpreted as the index of two-centered black hole 
solutions, carrying charges a.% = 71 and «2 = 272 for the first term, or a\ = 71+72 and 0,2 = 72 
for the second term. Indeed, for such two-centered solutions, the distance is fixed to 



r^2 = — - 
2 



1 (ax, 02) \Z(a\) + Z(a 2 ) 



Im [Z(a 1 )Z(a2)] 



(4.3) 



which is positive on the side c_ of the wall only; the unit vector ryilr\2 can be chosen arbitrarily 
on the unit two-sphere. Such configurations carry angular momentum J = \(a.\2 — 1)772/772 



(similar to the angular momentum carried by an electron in a magnetic monopole background), 
and therefore have |ai2 1 possible configurational states, with index (7(01,02) = (— l) ai2 oi2. 
Since the distance r\2 diverges in the vicinity of the wall, the internal degrees of freedom 
are decoupled from the configurational degrees of freedom, and the total index is the product 
g(a\, O2)f2 _ (oi)0~(o2), consistently with (4.2). 

The contribution on the second line of (4.2) is more interesting. Letting d = 7i2^ + (0, 1), the 
term in bracket is recognized as the index of the symmetric (or, when d < 0, antisymmetric) 
part of the tensor product Hi (g> Hi, where Hi is the space of quantum states accessible to one 
particle of charge 72 in the field of the particle of charge 71. In other words, the second line 
corresponds to a configuration of two identical centers of charge 72 orbiting around a center of 



charge 71, with Bose statistics when d > 0, or Fermi statistics when d < 0. More generally, (4.1 ) 



can be interpreted as the contribution from halos of particles of charge k~/2 orbiting on a fixed 
shell around a center of charge 71 + ^072, and obeying Bose or Fermi statistics, depending on 
the sign of (71 + /c 72, ^72) ^ + (^72)- 

In contrast, in terms of the rational DT invariants the semi-primitive wall-crossing formula 



(2.22) reads 



N 



so that, e.g. for 7 1— > 71 + 272, 



exp jT(-l)<<™> £(71,72) ^72, y) q l 



(4.4) 



,£=1 



Afi(i, 2) = 2 7l2 n + (i, 0) n+(o, 2) + (-lp 2 712 i)n + (o, 1) o) 



(712^(0,1))' 



where we combined f2(0, 2) and — |O(0, 1) in (4.2) into 0,(0,2). Unlike (4.2), the last term in 



(4.5) is of the form \d 2 , 



which would be the result if the two particles of charge (0, 1) were 



distinguishable and obeyed Boltzmann statistics. More generally, (4.4) can be interpreted as 



contributions of the same halo of particles with charge k^i described above, but now satisfying 
Boltzmann statistics. While (4.5) is hardly shorter than (|4.2|), the reader can easily convince 



him/herself of the power of this simplification by computing AQ(1,N) for higher N. 

The lesson to take from this discussion is that, rather than computing the variation of 
fi(7, t) across the wall W(7i,72), it is advantageous to compute instead the variation of the 



rational invariants 0(7, t) defined in (2.13), and apply the following recipe: treat the centers 
as distinguishable pointlike particles and compute their configurational index. Whenever m 
centers carry the same charge o, divide the configuration index by a Bolztmann-Gibbs factor 
1/m!. Finally, multiply the configurational index by the effective (rational) index Cl(a,t) carried 



by each of the centers. The result of this recipe is the formula (1.4), where g({a{\) is identified as 
the configurational index of the quantum mechanics of n distinguishable particles interacting via 
Coulomb and Lorentz type forces (which we discuss in detail in the next subsection). Although 
we motivated this recipe by inspecting the classical semi-primitive formula, it in fact holds in 
full generality and applies to the refined (or motivic) index as well, see [1] for more details. 



4-2. The phase space of multi- centered BPS black holes 

Let us now review some relevant properties of supersymmetric multi-centered black hole solutions 
in N = 2 supergravity (a similar analysis for multi-centered dyon solution in the low energy 
limit of N = 2 gauge theories can be found in [30|. I31j). Such solutions fall into the stationary 
metric ansatz 

ds 2 = _ e 2U {dt + A) 2 + e -2U df 2 (4g) 



where the scale function U, the Kaluza-Klein one- form A and the vector multiplet scalars 
z a , a = 1 . . . n v depend on the coordinate r on M 3 . 

For n centers located at r±, . . . , f n , carrying electromagnetic charges a±, . . . , a n G T with total 
charge 7 = a\ + ■ ■ ■ + a n , the values of the vector multiplet scalars t and of the scale factor U 
are obtained by solving (T2| . 



- 2 e~ u ^ Im \e-^Y(t(r))] = (3 + ^ 



i=l 



Til ' 



argZ(7,too), 



(4.7) 



where l^(t) = —e ' (X (t), F\(t)) is the symplectic section afforded by the special geometry of 
the vector multiplet moduli space, such that ^(7, t) = (7, y(too))- The constant vector (3 on the 
right-hand side of ( |4.7| ) is determined in terms of the asymptotic values of the moduli at infinity 

too by 



-21m 



(4.8) 



In particular, it follows from (4.7) that the scale factor U is given by evaluating the Bekenstein- 



Hawking entropy function ^(7) on the harmonic function appearing on the right-hand side of 

1 



(4.7) [32J 



-2U(f) 



7T 



n 

E 

8=1 



(4.9) 



Most importantly for our purposes, the locations r*j are subject to the condition of mechanical 
equilibrium under the Coulomb, Lorentz, Newton and scalar exchange forces (also known as 
integrability equations) [12] 



E 

3=1 



a 



X) 



(4.10) 



where r 



'j 



{ai,a.j), and the real constants 



depend on the the asymptotic values of the moduli. Since 



(4.11) 

: arg Z(7, too), the constants q 
satisfy Y27=i Ci = ®- ^ ne conditions (4.10) guarantee the existence of a Kaluza-Klein connection 
A such that the above configuration is a supersymmetric solution of the equations of motion. 
In order for the solution to be physical however, one must also require that the scale factor be 
everywhere positive 



S 



> 0. 



V r G 



(4.12) 



i=l 



where r*j is the location of the i-th center. For the configurations relevant to the wall-crossing 
problem, this condition appears to be automatically satisfied. 

Now, let Ai n ({ctij}; {q}) be the space of of solutions {n, . . . f n } to the equilibrium conditions 
(4.10), modulo overall translations of the centers. Ai n is a (possibly disconnected) 2n — 2- 
dimensional submanifold of M 3n_3 \A, where A is the locus in M 3 ™ -3 where two or more of the 
centers r*j coincide. M 3 ™ _3 \A is equipped with the closed two-form 



U) 



1 „ e abc dr* A drf, rf, 
4 2^"^ ^..13 



(4.13) 
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Figure 2. Phase structure of the moduli space M n of 3-centered solutions as a function of Cj, 
for fixed charges such that a%2 > 0, «23 > 0, a.\% > 0, a\i < 023- The shaded area represent the 
values of Cj in the two-dimensional section c\ + C2 + C3 = which are spanned as the location of 
the 3rd center is varied, keeping the centers 1 and 2 fixed. Conversely, if the values of Cj is fixed, 
the range of distances between the centers 1 and 2 can be read off by intersecting the shaded 
area with a radial line which joins Cj to the origin. Thus, for this choice of charges, 3-centered 
solutions only exist in the region c\ > 0, C3 < 0. Inside this region, the range of 7*12 is bounded 
from below and from above, except on the wall of marginal stability C2 = 0. The boundaries 
of the shaded region correspond to collinear solutions whose order is indicated. As the wall is 
crossed, the topology of the collinear solutions changes from (321), (132) to (321), (312). 

For generic values of Cj, the restriction of uj to M n is non-degenerate and endows M n with a 
symplectic structure [16] . Moreover, the symplectic form cj is invariant under 50(3) rotations. 
The moment map associated to infinitesimal rotations is the angular momentum 



Away from walls of marginal stability, the distances are bounded from above. If it is possible 
to order the ctj's such that (ati,aj) > whenever i < j, as it is the case when the a^s lie 
in a two-dimensional cone T, the distances rjj are also bounded from below by a non-zero 
^min > 0, and the space M. n is therefore compact. E.g. for two centers, M.2 = S 2 equipped with 
oj = ^a\2sm9 d8d(j) for sign(ai2) = sign(ci), and zero otherwise. A representative of the phase 
structure of M.3 is illustrated in Figure [2] 

4-3. Equivariant index and localization 

The symplectic space Ai n = M. n ({aij}; {ci}) defines the classical phase space of the 
configurational degrees of freedom of n-centered BPS black hole solutions. Since such 
configurations are stationary, the Hamiltonian vanishes and all points in Ai n are degenerate 
in energy. Quantum mechanically, the Hilbert space consists of sections of % = S <g> £, where 




(4.14) 



S = S+ © S- is the total spin bundlfl over Ai n and C is a complex line bundle over Ai n with 
first Chern class uj. BPS states correspond to zero-modes of the Dirac operator D on Ai n , and 
are analogous to states in the lowest Landau level for an electron immersed in a magnetic flux 
uj. The Dirac operator decomposes as D = D + + D_ where D + maps S + © C to 5_ © £ and 
vice-versa. The action of 5*0(3) on .M n lifts to an action of SU(2) on S± © £, and the refined 
index is then 

v) = Tr Ke rD + (-y) 2j3 + Tr Ker D_(-y) 2j3 (4.15) 

where J3 is the operator representing the rotations along the z axis. Assuming that KerZ?_ = 0, 
it follows that the refined index is equal to the equivariant index of the Dirac operator D, 

g Ie{ ({ ai },y) = Tr KcrD+ (-y) 2J * - TV K erD_(-y) 2J3 (4.16) 

The assumption that KerD_ = can be proven when A4 n is Kahler (which is the case when 
n = 2,3). We do not know how to prove it in general, but it is supported by the fact that it leads 
to results in agreement with the KS and JS formulae. In <|6j we speculate that this assumption 
may be unnecessary if one were to compute the jump of the protected spin character in the 
context of M = 2 SYM theories. 

Now, by the Atiyah-Bott Lefschetz fixed point formula |33 | 134 } [35 | I36j. the equivariant index 
localizes to the fixed points of the action of J3 on M. n - Clearly, those correspond to solutions 
where all centers lie along the z axis, and satisfy the one-dimensional equilibrium conditions 



j=i ^ Zti i=i 

where the last equation fixes the translational zero- mode. For any a 6 S n , where S n denotes 



the set of permutations of {l,...,n}, we denote by C(a) the set of solutions to (4.17) such 
that z a u\ < Zcr{j) if i < j- The set C(a) corresponds to the subset of the critical points of the 
'superpotential' 

W(X, {zi}) — - aij sign(zj - Zi) \n\zj - z%\ - ^(cj - \/n)zi (4.18) 

i<j i 

which are ordered according to the permutation a. In the vicinity of a fixed point p G C(a), the 
angular momentum J3 and the symplectic form oj take the form 

^3 = 2 ^2 a <r{i)<r<J) ~ T^'iW ( X i X 3 + ViVj) H > U = \ M ij(p) dx i A d Vj H > ( 4 - 19 ) 

i<j 

where Mij is the Hessian matrix of W(X,{zi}) with respect to zx,...z n , and (xi,yi) are 
coordinates in the plane transverse to the z-axis at the center i, subject to the condition 
Si x i = S Vi = 0- Except for an overall translational zero-mode, the matrix My is non- 
degenerate, and the critical points are isolated, so C(a) is a finite set (possibly empty). 
The Lefschetz fixed point formula of yields an explicit formula for the refined index 



(■_l)Ei<iOij+n-l 

(y - 1/y) 1 



4 We assume that A4 n is spin, see [TT] for a discussion of this issue. 



where s(a) counts (with sign) the number of solutions to (4.17) ordered according to the 
permutation a, 

s(a) = - signdetM, (4.21) 

pec(a) 

where M is the Hessian of W a,s ci function of the ti ~\~ 1 variables A, . . , Zn-> evaluated at the 
given solution of (4.17). The factor (y — l/y) n in (4.20) originates form the equivariant A- 



genus in the Lefschetz fixed point formula. It is convenient to let the sum in (4.20) run over all 
permutations and set s(a) = when there are no solutions to (4.17) in the order specified by a. 



For reasons that will become clear in we refer to (4.20) as the 'Coulomb branch wall-crossing 
formula'. 



The formula (4.20) is fully explicit, yet it depends on our ability to find solutions of the 



one-dimensional problem (4.17). While this can be done numerically (approximate solutions are 
sufficient since the answer depends only on the order a and the sign of detM), it would be 
useful to have a general criterium for determining when solutions exist, and if so to compute 
their Morse index. The answer to these questions is suggested by an alternative approach based 
on quivers (see ^6|. 



At this point, we can check whether (4.20) agrees with the answer of the KS or JS formulae. 
For n = 2, we find two fixed points with permutations cr (12) = 12 and 21, leading to 



g(ai,a 2 ;y) 



, sinh(z/ai2) 
sinh v 



(4.22) 



in agreement with (2.19). For n = 3 (and for the same choice of a.% as above (2.20)), we find 4 
possible orderings, with the following value of s(a), 



{1,2, 3;+}, {2, 1,3;-}, {3, 1,2;-}, {3, 2,1;+} 



(4.23) 



leading to 

SrefCai.aa.ag.y) = (-l)^+^+a 13 (y _ y -iy2 

CI12+ai3+a23 _ „,«13+ a 23— a 12 _ ,,012— «23~ "13 



x \y 



v 



v 



+ y 



-Ct\2— «13— "23 



(4.24) 



in agreement with the result (2.20). Further checks for n > 3 can be found in pQ. It is an open 



problem to show by combinatorial means that the result (4.20) agrees with the KS formula in 
general. 



5. Wall-crossing from Abelian quivers 

Finally, we discuss an alternative computation of the combinatorial factors g re { ({ai}, y) based 
on Reineke's formula for invariants of quivers without closed loop. 

The physical idea is that the supersymmetric quantum mechanics of multi-centered black 
holes admits two branches: the Coulomb branch, where the centers are far separated and well 
described by the supergravity solution described in §4.2[ and the Higgs branch, where the centers 
are very close to each other and are better represented as D-branes, with open strings stretched 
between them. At large string coupling, the wave function is mainly supported on the Coulomb 
branch, while at small string coupling it is mainly supported on the Higgs branch [JS]. However, 
the BPS index Tr'(— 1)^ is independent of the string coupling (and other hypermultiplet fields), 
and should be computable in both. It is less clear that the same should be true of the refined 
index Tr'(— y) 2Ji , since this quantity is in general non-protected in string theory [28], but for 



what concerns the jump of the refined index across walls of marginal stability, we shall find 
strong evidence that this is the case. 

On the Higgs branch, the D-brane system is described at low energy by a supersymmetric 
quiver quantum mechanics, with gauge group Y]a=i ^C^i)j wriere Ni is the number of coinciding 
D-branes at point i, and (oj, a,- fields in the bifundamental representation (iVj, Nj) when (oj, otj) 
is positive, or —{on, (Xj) fields in the bifundamental representation (iVj, Nj) in the opposite case. 
Thanks to the Bose-Fermi/Boltzmann correspondence described in £ |4.1[ we can treat all centers 
as distinguishable and assume that N = 1 for any i, provided that we attach an effective rational 
index Q ie f(ai,y) to each center, and perturb the charge vectors so that none of them coincide. 
Moreover, since the on can be ordered such that (oj,Oj) > whenever i < j, the quiver admits 
no closed loop. 

For arbitrary quivers without loops, Reineke has computed the Poincare polynomial of the 
moduli space of quiver representations [37]. In the special case of Abelian quivers (N = 1), 
Reineke 's formula gives 



({«t},y) 



-y) 



{y-l/yY- 1 4-* 

partitions 



\ " f_iy-i y 2 '£a<b'Ej < i a ji m i a)m j' ) 



(5.1) 



where the sum runs over all ordered partitions of 7 = a\ + ■ ■ ■ + a n into s vectors (3^ (1 < a < s, 
1 < s < n) such that 



(i) = ^mf'ai with mf> G {0, 1}, Y. a P {a) = 7 

(ii) (ELi P {a \l) > v b with l<b<s-l 



We refer to (5.1) as the 'Higgs branch wall-crossing formula'. 



Let us illustrate how this formula works by computing g re { ({cti}, y) for n = 2,3. For n = 2 
case with 012 > 0, there are two possible ordered partitions satisfying the conditions stated 
above: 

{ai + a 2 }, {ai,a 2 }. (5.2) 
The first term contributes y 2 " 12 and the second term contributes —1 to the sum. In total, 



9 ref(ai,a 2 ,y) = (-y) 1 ^ 12 (y 2 - l)' 1 (y 2a ^ - 1) = (-l)^ 1 V - 



a l2 _ y~ ai2 



y-y 



(5.3) 



in agreement with (2.19). For n = 3, assuming the same conditions on a, as in (2.20), we find 



6 possible ordered partitions 

{01+02 + 03}, {01,02 + 03}, {01 + 02,03}, {01 + 03,02}, {01,02,03}, {01,03,02}. 

(5.4) 

The second and the last contribution cancel, leaving 



ffref («1, «2, « 3 , V) =(-l)^+«i3+a 23 {y _ y -l) 



-l\-2 



yCt 1 2+a 1 3+a 2 3 _ yOl2-a23-"13 _ yOl3+«23-"12 _|_ y~ "12" Ol3-«23) 



(5.5) 



in agreement with (2.20). 



6. Conclusion and open problems 

In this survey, we have described four apparently different wall-crossing formulae, the KS 



formula (2.9), the JS formula (3.6), the 'Coulomb branch' formula (4.20) and the 'Higgs branch 



formula (5.1). In [TJ we have checked that these formulae agree among each other for n < 5. 



Unfortunately, we do not have a mathematical proof of their equivalence at this point. The 
equivalence of the KS and JS formula appears to be on solid ground [TUJ, [T| . The Reineke formula 
(or rather, its specialization at y = 1) is known to follow from the JS formula [3H] (although I do 
not know of an explicit combinatorial proof). The underlying mathematical structure of these 
formulae involves Ringel-Hall algebras (and generalizations thereof [39]), which might provide a 
realization of the long sought-after algebra of BPS states [101 HI] . 



The Coulomb branch formula (4.20) appears to be new, and it would be desirable to derive 



it e.g. from the Higgs branch formula (5.1 ). Indeed, the equivalence between these two formulae 
could be viewed as a toy model of open/closed string duality. The comparison of the Higgs and 
Coulomb branch formulae suggests that the permutations a for which s(a) does not vanish are 



those whose maximal increasing subsequence ^] satisfies the condition ii) below (5.1), while none 



of its non-maximal increasing subsequence does, in which case s(a) = (— l)#i^°T i + 1 ) <cr (*)}. It 
would be interesting to derive these conditions by applying Morse theory to the 'superpotentiaP 



W in ( [4.18D . 

While the Coulomb and Higgs branch formulae have been derived in the context of J\f = 2 
supergravity, the fact that they appear to agree with the KS and JS formulae suggests that they 
should work just as well in the context of N = 2 SYM theories. For such field theories, unlike 
in supergravity, the refined index (or rather a variant thereof, known as the protected spin 
character, and constructed out of the spatial rotation generator J3 and a SU(2)r symmetry 
generator I3 |28j ) receives contributions from short multiplets only. It would be desirable to 
derive the jump of the protected spin character by quantizing multi-centered dyonic solutions 
of the low energy Abelian gauge theory, along the lines of |3Q|. 131] , and see if it agrees with the 
Coulomb formula (presumably this jump will coincide with the equivariant index of the Dirac 
operator D, without the need to assume that KerD_ =0). It would also be interesting to clarify 
the relation with other realizations of BPS dyons in N = 2 gauge theories based on string webs 
or non- Abelian monopoles [4"2 l |4"3 1 1311 1351136] . 

More generally, localization techniques appear to be a powerful way of quantizing multi- 
centered black hole solutions at fixed values of the moduli, not only in the vicinity of walls of 
marginal stability. In contrast to the situation studied here, the phase space Ai n is in general 
no longer be compact, due to the presence of scaling solutions [T31 [371 HH1 HH HH]. We refer the 
reader to [TT] for an application of localization techniques to this problem. 
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